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Abstract 

It has been recently conjectured by Selem and Wilczek [1] the existence of a 
ss — [ud] meson due to strong correlations between the two hght antiquarks. 
We make a detailed study of this system within a dynamical quark model 
which has proven to be successful in reproducing the most important fea- 
tures of low-energy hadron phenomenology. Our results, obtained within a 
parameter-free calculation, show that the antidiquark component of the ssud 
system indeed entails the stronger attraction, and drives its energy much lower 
than the Nr. threshold, but still above the K*~ or K*^ K~ thresholds. 
We have also studied the ccud and bbud systems. Exotic mesons are only ex- 
pected to exist in the limit of large mass for the two-quark subsystem, bbud, 
since the calculated mass is below the B° B*~ or B*'^ thresholds. 
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It has been recently re-emphasized [1-3] the potential importance of strong diquark (an- 
tidiquark) correlations in hadronic physics [4,5]. Theoretically the idea of diquark (antidi- 
quark) correlations inside hadrons is a consequence of color cancellation. The disturbance 
produced by the color charges of two quarks in empty space can be halved by bringing 
them together into a single 3 representation of the color SU{3) group. If this is joined 
with the more favorable spin-singlet state and Fermi statistics, the quarks must be in the 
antisymmetric 3 representation of flavor iS'?7(3). Besides, one should expect that any effect 
disrupting the correlations will induce a repulsive force. Among such effects we may quote 
the presence of an additional diquark or a spectator quark. Such ideas suggest that the eas- 
iest way of constructing low-energy exotics could be based on strongly correlated diquarks 
(antidiquarks) as building-blocks. 

However, these arguments are rather qualitative, and merely based on the group theo- 
retical structure of QCD. To study quantitatively whether or not such strong correlations 
between the light quarks (antiquarks) are indeed present, one needs QCD-based dynamical 
studies such as lattice QCD, although this eventually must be checked by experiments. In 
view of present status of the lattice QCD simulations [6], it is still meaningful to use phe- 
nomenological models which contain the main features of QCD, once the model parameters 
are calibrated and constrained by as many observables as possible. Here, we perform such 
a consistent, parameter- free dynamical calculation for the study of the correlations between 
light antiquarks, and investigate the possible existence of the exotic meson, ss — [ud\ , which 
has been recently conjectured by Selem and Wilczek [1]. For completeness, we have also 
analyzed the ccud and bbud systems. 

The present study has been done within the framework of a constituent quark model 
which has been successfully applied to study the baryon spectra and the baryon-baryon 
interaction [7]. This model has been generahzed to include also strange (s), charm (c), 
and beauty (6) flavors, and it has also been shown to give a reasonable description of the 
meson spectra [8]. The description of experimental data gets improved when four-quark 
{qqqq) components are also considered [9,10]. The model parameters have been strongly 
constrained by the study of different hadron observables, what represents an advancing 
feature compared to studies based on models designed adhoc for a particular problem. 

The model is based on the assumption that the u, d and s constituent quarks acquire their 
masses due to the spontaneous breaking of the original SU (3) l ® SU chiral symmetry 
at some momentum scale, which is one of the most important nonperturbative phenomena 
for low energy hadron structure. In this domain of momenta quarks are quasi-particles with 
constituent masses interacting through scalar (sigmas, OSE) and pseudoscalar (pions, OPE; 
kaons, OKE; and etas, OEE) boson-exchange potentials. Note that for the case of heavy 
quarks, c and b, boson-exchange potentials are not present in the model [10], since chiral 
symmetry is badly broken already at the level of the current quark masses. Beyond the 
chiral symmetry breaking scale one expects the dynamics being governed by QCD perturba- 
tive effects. They are taken into account through the one-gluon-exchange (OGE) potential, 
a standard color Fermi-Breit interaction. Finally, any model imitating QCD should incor- 
porate conflnement (CON). Lattice calculations in the quenched approximation for heavy 
quarks show that the conflning interaction is linearly dependent on the interquark distance. 
The presence of sea quarks, apart from valence quarks (unquenched approximation), sug- 
gests a screening effect on the potential when increasing the interquark distance. Creation 
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of light-quark pairs out of vacuum in between the quarks becomes energetically preferable, 
resulting in a complete screening of quark color charges at large distances. String breaking 
has been definitively confirmed through lattice calculations [11] in coincidence with the quite 
rapid crossover from a linear rising to a fiat potential well estabhshed in SU (2) Yang-Mills 
theories [12]. Explicit expressions for the interaction potentials derived from the nonrela- 
tivistic reduction of the Lagrangian density in the static approximation, and a more detailed 
discussion of the model can be found in Ref. [8]. 

For the description of the most general QQud {Q = s, c, or b) system we introduce the 
Jacobi coordinates, 

f = fi-f2, y = r3-u, 
z = , 

mi + 7722 7723 + 7714 

where 1 and 2 (3 and 4) stand for quarks (antiquarks). The ground state energy of the 
four-body problem can be estimated by a variational method using a trial wave function 
that includes all possible color-flavor-spin components relevant to a given configuration. For 
each component, | 0i >, such a basis wave function will be a tensor product of color (cj), 
fiavor (/j), spin [xi) ^'^^ spatial [Ri) parts, 

I 0i >=| Ci(1234) > ® I /i(1234) > ® I Xi(1234) > ® | i?i(1234) > . (2) 

The most general spatial wave function can be expressed as a combination of six scalar 
quantities, 

I i?,(1234) >= R,{x\y\z\x-y, x-z,y-z}. (3) 

The variational spatial wave function is taken to be a linear combination of generalized 
Gaussians, 

I it:,(i234) >= f: f^^E^i^ = x: (5^'^^^'^'-^'y'-'^'''-^'^-y-<^'^-^^^^^^ (4) 

where n is the number of terms to expand the spatial wave function of each color-fiavor-spin 
component, and a:p ^ b'f\ /P'* are the variational parameters. 

With respect to the color wave function, | Ci(1234) >, one can couple the two quarks 
(1, 2) and the two antiquarks (3, 4) to a color singlet state in different ways: 

I li3, 124 > , I 813,824 > ; (5) 
I Ii4,i23 > , I 814,823 > ; (6) 

I 3i2 , 334 > , I 612,634 > . (7) 

The couplings in Eqs. (5) and (6) are convenient for asymptotic meson-meson channels (or 
meson-meson molecules) while those in Eq. (7) are more appropriate for tetraquark bound 
states. With our choice of the Jacobi coordinates the color basis in Eq. (7) is more suitable 
to deal with the Pauli principle in an easier way. 
The spin part can be written as. 
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X.(1234) >= [{12)s,,mss,]s, 



(8) 



where the spin of the two quarks is coupled to 5*12 and that of the antiquarks to 5*34. 

Concerning the flavor part, | /j(1234) >, since the heavy quarks (those with flavor s, c or 
b) have isospin zero, they do not contribute to the total isospin. Therefore one can classify 
the tetraquark wave function by the isospin of the light quarks / = 0, 1. Taking into account 
all degrees of freedom, the Pauli principle must be satisfied for each subsystem of identical 
quarks (antiquarks). It restricts the quantum numbers of the basis states, that justifies to 
use the [{QQ){ud)] coupling. 

Using the wave functions described above, we search for a variational solution for the 
Hamiltonian. The color, fiavor and spin parts are integrated out and the coefficients of 
the spatial wave function are obtained by solving the system of linear equations, 

n 

E E [{^ I H \B!f ) - E {lf^\B!f %,,] = for aU k, t', (9) 

i j=l 

once the eigenvalues E{a\-'\ b\-'\ //"'■') are obtained by a minimization procedure with re- 
spect to the variational parameters. The stable tetraquark states are identified by comparing 
the obtained eigenvalues with the corresponding physical thresholds. If they are above the 
threshold they would be very broad objects, very hard to detect experimentally. 

In a realistic model tetraquarks will not overpopulate the meson spectra, in fact they 
may complement two-quark components and, indeed, they seem to be necessary in order 
to understand the rich meson phenomenology [9,10]. This is due, on the one hand, to 
the constituent mass of the quarks, and on the other one, to the finite spectra generated 
by screened confining potentials [13]. Only positive parity tetraquark states, those that 
do not need internal orbital angular momentum between the constituents, may appear in 
the low-energy region of the meson spectra and they could mix with qq states with the 
same quantum numbers. Negative parity four-quark states need a unit of orbital angular 
momentum what means an average excitation energy of 800—900 MeV [14]. These ideas have 
been recently used to explain the abnormal number of low-energy scalar-isoscalar mesons [9] 
and also the unexpected low masses of positive parity (O"*" and 1"*") open-charm mesons [10]. 
They are perfect examples of the way how the enlargement of the Fock space may help in the 
understanding of meson phenomenology. As explained in these works, only those states with 
exotic quantum numbers may appear as pure four-quark resonances on the meson spectra. 
Unfortunately, the present uncertainties on the experimental data concerning exotic channels 
prevents, for the moment, to extract a definitive conclusion about its existence [15]. 

Let us concentrate on the particular meson state, ss — [ud], conjectured in Ref. [1]. It 
has the property of the [ud\ subsystem being an antidiquark state, which means the two 
antiquarks are in a color (Sg), fiavor {1^^^] = 0); spin (^jj^jj = 0) antisymmetric state. 
It requires a completely symmetric radial wave function for the two antiquarks to satisfy 
the Fermi statistics. As antidiquark component [ud] should be in a relative S"— wave, one 
can neglect the crossing terms in the trial radial wave function, those depending on the 
scalar product of different Jacobi coordinates in Eq. (4). In order to obtain a color singlet 
wave function, the two s quarks must be in a color antisymmetric, 3c, state. Being flavor 
symmetric, the corresponding spin wave function may be either in (i) an antisymmetric. 
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Sss = 0, state that would require the anti-natural radial antisymmetric wave function to 
describe the ground state of the system, or in (ii) a symmetric spin state, Sgs = 1, that 
would combine with a natural symmetric radial wave function. Therefore, the conjectured 
meson with the presence of the antidiquark would be described byaJ'^ = l+(L = 0,5'=l) 
state with isospin 7 = 0. We can summarize the quantum numbers of the antidiquark 
component of the ssud system in the following way, 

I [3e, S = 0,1 = QUA [3c, S=l,I = OUiS=l,I = 0)). (10) 

A full calculation of the J'^ = 1"*" (L = 0, 5" = 1) state with strangeness —2 would require 
also to consider other vectors in the Hilbert space. In particular, the same state could also be 
constructed from a different vector, where the two antiquarks would not be an antidiquark 
state while it still has a completely symmetric radial wave function, 

I [6e, S=1,I = 0]^j, I [6„ S = Q,I = OU {S=1,I = 0)). (11) 

This vector, which will be referred to as the nondiquark component of the ssud system, 
should be considered in the calculation of the four-quark state without requiring the antidi- 
quark configuration, whereas, it will not be included if only the antidiquark configuration is 
imposed. For the full calculation, both the antidiquark and nondiquark configurations will 
be included with the corresponding configuration mixing. 

In Table I we present our results for the antidiquark configuration, nondiquark configu- 
ration, and full calculation for the ssud system. The same calculation has been repeated for 
the ccud and bbud systems, and the results are presented in Tables II and III, respectively. 

The first important conclusion that can be extracted from the results of Tables I, II and 
III is that the energy of the antidiquark configuration is always the lowest. It is interesting 
to note how the pseudoscalar force acting between the light quarks is responsible for that, 
since the results for the antidiquark and nondiquark configurations are almost degenerate if 
only the confinement and one-gluon exchange are retained. The reason for this stems on the 
different symmetry for both components in color-spin and fiavor-spin spaces. While both 
are symmetric in color-spin space, the antidiquark (nondiquark) component is symmetric 
(antisymmetric) in flavor-spin space. Therefore, if strong diquark correlations were dictated 
by QCD for light quarks, the dynamical explanation could not rely on the simple one- 
gluon exchange dynamics, but it would need meson-exchange forces between the constituent 
quarks. The similar effect has been also observed in the case of baryon spectra, where 
pseudoscalar meson exchanges between the constituent quarks are able to revert the relative 
position in the energy spectra of the nucleon Roper resonance, with a dominant flavor- 
spin symmetric wave function, with respect to negative parity states, with a flavor-spin 
antisymmetric wave function [16]. It is also interesting to notice that the antidiquark and 
nondiquark components are not exactly degenerate when only the confining interaction is 
taken into account. This can be easily understood by looking at Table IV, where we present 
the contribution of the interaction between QQ, V12, ud, V34, and Qn {n = u,d), Vis, for the 
QQud system as a function of the mass of Q, mq. The minimization procedure modifies the 
variational parameter a^p in Eq. (4) for the coordinate due to the smaller size of the QQ 
subsystem when the mass mq increases. As a consequence it gives a smaller contribution 
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to the energy of the system. In other words, the dependence on the mass of the quark is 
introduced into the calculation through the variational parameters. 

As predicted by Selem and Wilczek [1], the mixing between the antidiquark and nondi- 
quark components of the wave function diminishes when increasing the mass of the heavy 
quarks (see Table V) , in such a way that for the b quark case the nondiquark component gives 
almost no contribution to the ground state energy of the system. This effect, interpreted 
as the less capacity of the spin of a heavy quark to disrupt the correlation of the diquark 
(antidiquark), comes from the l/(mjmj) {rriij: constituent quark masses) dependence of 
the one-gluon exchange interaction, which is responsible for the mixing between these com- 
ponents. Therefore, the mixing decreases with increasing the mass of the heavy quark. A 
similar evidence, the less capacity of the spin of the heavy quarks to disrupt the system, 
has also been observed in the spin-orbit splitting of the A— , A+— and A;,— hypernuclei [17], 
where s, c and b quarks exclusively carry the total spin of the A, A+ and A^, respectively, 
and u and d quarks are coupled to a isospin zero and spin zero diquark in each baryon. 

Regarding the possibility of observing these systems, the results obtained are always far 
above their corresponding lowest two-meson thresholds, as indicated in Tables I, II, and III, 
being the only exception the bbud system. 

Experimentally the possibility to detect a QQud meson relies on two different aspects. 
First of all the rate of production of QQud, and second the existence of decay modes that 
can provide a unique signature. For the production at hadronic or e+e~ colliders one needs 
to produce two pairs of charm or bottom quarks. These pairs should be close spatially and 
the quarks within each quark-antiquark pair should have small relative momenta in order 
to combine in a two-quark, cc or bb, state. Finally, these two-quark states should pick up an 
antidiquark [ud\ to form the desired QQ — [ud] system. The production rates for the case of 
charm quarks have been estimated in Refs. [18] and [19]. The signal of the strong antidiquark 
correlation would come from decay channels preferring to keep the antidiquark structure. 
So, instead decaying by splitting into a two-meson system, it would proceed through a two 
baryon system as it would be NScc for the charm case and NEn for the bottom case. The 
absence of a dynamical enhancement of the antidiquark component would open the decay 
into two mesons. For the ssud system, if the antidiquark component is strong enough so as 
to force it to decay into a two baryon system via {ss — [ud]) NE, one can expect a natural 
decay in an S*— wave, which needs a J'^ = 1~ state for the system. As the energy excitation 
for a unit of orbital angular momentum costs about 800—900 MeV [14], this would make the 
system to be in the continuum above the NE threshold, therefore being broad and difficult to 
detect. The dynamical enhancement of the diquark component is one of the possible reasons 
to explain the decay of the A(3/2~)(1520) 45% of the time to NK channels in order to retain 
the diquark structure of the A inside the N. Finally, let us remark that we only support the 
possible existence of QQud states for Q = b and with more uncertainty for Q = c, but never 
for Q = s. In the bottom case the predicted state should be very narrow and easy to observe 
(if produced) since it is far below the two-meson and two-baryon physical thresholds. 

To summarize, we have made a dynamical, parameter-free, calculation for the QQud 
system {Q = s, c, b) within a realistic constituent quark model. We have found that the 
antidiquark configuration of these systems always gives a lower energy than the nondiquark 
configuration. The mixing between the antidiquark and nondiquark states is due to the 
one-gluon exchange potential, and because of its l/{mimj) dependence, it decreases when 
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increasing the heavy quark mass. As a consequence one can expect that the conjectured 
mesons QQ—[ud\ could be stable for Q = b, but we do not find any reason why these systems 
should be bound for Q = s. Moreover, for light quarks, Q = s, there is no dynamical reason 
why the antidiquark component should be favored compared to the nondiquark one (within 
an uncorrelated quark model). If the existence of any such systems with two light quarks, 
Q = s, and therefore with a strongly correlated light antidiquark would be confirmed via 
the postulated baryon-antibaryon final channel, it would mean that a dynamical mechanism 
responsible for the correlations has not been considered in our simple realizations of models 
for QCD (based on uncorrelated quarks), and the question for the existence of exotic systems, 
such as the pentaquark, should be addressed in a corresponding manner. On the other hand, 
the present results may imply that the pentaquark with a heavy b quark, [ud] [ud]b with a 
negative parity, may have a chance to be stable, if the predicted repulsion between diquarks 
is not strong enough to destroy the system. This has a merit of further study within the 
same model. 



I. ACKNOWLEDGEMENTS 

This work has been partially funded by Ministerio de Ciencia y Tecnologia under Contract 
No. FPA2004-05616, by Junta de Castilla y Leon under Contract No. SA-104/04, and by 
Generalitat Valenciana under Contract No. GV05/276. 



7 



REFERENCES 



[1] A. Selem and F. Wilczek, hep-ph/0602128. 
[2] R.L. Jaffe, Phys. Rep. 409, 1 (2005). 

[3] R. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003). 
[4] H.J. Lipkin, Phys. Lett. B 35, 534 (1971). 
[5] M. Karhncr and H.J. Lipkin, hcp-ph/0601193. 

[6] E.G. Lasscock et al, Nucl. Phys. B, Proc. Suppl. 153, 348 (2006). H. Suganuma, K. 

Tsumura, N. Ishii, and F. Okiharu, PoS LAT2005, 070 (2005), hep-lat/0509121. 
[7] A. Valcarce, H. Garcilazo, F. Fernandez, and P. Gonzalez, Rep. Prog. Phys. 68, 965 

(2005). 

[8] J. Vijande, F. Fernandez, and A. Valcarce, J. Phys. G 31, 481 (2005). 
[9] J. Vijande, A. Valcarce, F. Fernandez, and B. Silvestre-Brac, Phys. Rev. D 72, 034025 
(2005). 

[10] J. Vijande, F. Fernandez, and A. Valcarce, Phys. Rev. D 73, 034002 (2006). 

[11] SESAM Collaboration, G.S. Bah, H. Neff, T. Diissel, T. Lippert, and K. Schilhng, Phys. 

Rev. D 71, 114513 (2005). 
[12] P.W. Stephenson, Nucl. Phys. B 550, 427 (1999). P. de Forcrand and O. Philipsen, 

Phys. Lett. B 475, 280 (2000). 
[13] J. Vijande, P. Gonzalez, H. Garcilazo, and A. Valcarce Phys. Rev. D 69, 074019 (2004). 
[14] N. Barnea, J. Vijande, and A. Valcarce, Phys. Rev. D 73, 054004 (2006). 
[15] A.P. Szczepaniak, M. Swat, A.R. Dzierba, and S. Teige, Phys. Rev. Lett. 91, 092002 

(2003); A.R. Dzierba et al, hep-ex/0510068. 
[16] H. Garcilazo, A. Valcarce, and F. Fernandez, Phys. Rev. C 64, 058201 (2001). 
[17] K. Tsushima, K. Saito, J. Haidenbauer, and A.W. Thomas, Nucl. Phys. A630, 691 

(1998). K. Tsushima, and F.C. Khanna, Phys. Rev. C 67, 015211 (2003); Prog. Theor. 

Phys. Suppl. 149, 160 (2003); J. Phys. G 30, 1765 (2004). 
[18] B.A. Gelman and S. Nussinov, Phys. Lett. B 551, 296 (2003). 
[19] M.A. Moinester, Z. Phys. A 355, 349 (1996). 



8 



TABLES 



TABLE I. Calculated energies and physical thresholds, in MeV, for the ssud system. "Full" 
stands for the results calculated including botli configurations, antidiquark and nondiquark. Inside 
the brackets is the percentage of the antidiquark component in the full calculation. 





Antidiquark 


Nondiquark 


Full 


VcON + VoGE + VoPE + VoSE + VqEE 


1705 


1974 


1696 (97.59 %) 


VcON + VoGE + Vqpe + Vqse 


1644 


1965 




VcON + Vqge + VoPE 


1843 


2105 




VcoN + Vqge 


2092 


2083 




VCON 


2520 


2479 




N :l threshold 


2257 






A " K or K* K threshold 






1386—1390 


TABLE 11. Same as 


Table I for the ccud system. 






Antidiquark 


Nondiquark 


Full 


VcON + VoGE + VoPE + Vqse + Vqee 


3929 


4207 


3927 (99.48 %) 


VcoN + Vqge + Vqpe + Vqse 


3858 


4210 




VcoN + Vqge + Vqpe 


3906 


4229 




VcoN + Vqge 


4169 


4197 




VcON 


4631 


4644 




N threshold 


4460 






D+ D*" or D*^ t hroshold 






3875-3876 
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TABLE III. Same as Table I for the bbud system. 





Antidiquark 


Nondiquark 


Full 


VcON + VoGE + Vqpe + VoSE + VqeE 


10426 


10797 


10426 (99.95 %) 


TT" iT/ iT/" iT/ 

VcON + VOGE + VOPE + VqSE 


lUooo 


10801 




VcoN + Voge + Vqpe 


10403 


10822 




VcON + VoGE 


10673 


10787 




VcON 


11154 


11234 




N Sfeb threshold 

B*^B- or WB*' threshold 






10604 



TABLE IV. Expectation value, in MeV, of different contributions of the confining interaction 
for the different components of the QQud system and for two different values of the mass of the 
quark in the two-quark subsystem. 





nT-Q = 


555 MeV 


= 


5100 MeV 




Antidiquark 


Nondiquark 


Antidiquark 


Nondiquark 


< Vu > 


+458 


-263 


+241 


-176 


< V34 > 


+527 


-284 


+513 


-264 


< ^13 > 


+259 


+633 


+225 


+536 



TABLE V. Probability, in %, of the antidiquark component, QQ — [ud], as a function of the 
mass of the quark, in MeV, in the two-quark subsystem. We also give the mass of the QQud system 
in MeV. 





^QQud 


p(Qg - [ud]) 


313 


1431 


96.09 


555 


1696 


97.59 


755 


1980 


98.18 


1255 


2815 


99.11 


1555 


3352 


99.37 
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